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ABSTRACT 

Here we review some of the main issues related to 
multi-wavelength source identification and characteriza- 
tion, with particular emphasis on the field of X-ray sur- 
veys carried out over the last years. This complex and 
time-consuming process is going to represent one of the 
main difficulties over the coming years, when signifi- 
cantly larger surveys, both in area and depth, will be car- 
ried out with the new generations of space- and ground- 
based facilities like e.g. eROSITA, WISE, VISTA, Pan- 
STARRS, and LSST. The Virtual Observatory can offer 
a reliable way to approach to a new concept of data han- 
dling and multi-wavelength source characterization, pro- 
vided that uniform and rigorous data analyses and exten- 
sive quality checks are performed. 
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1. INTRODUCTION 

The increasing number of large surveys carried out over 
the last decade has revolutionized the fields of observa- 
tional cosmology and astrophysics but, at the same time, 
it has posed several problems. First of all, despite former 
surveys, which were primarily focused on censing and 
determining the properties of a given source class in a 
given band, recent surveys have been truly characterized 
by an observational multi-wavelength approach, whose 
final goal is to estimate the physical source properties 
over the entire electromagnetic spectrum and investigate 
how different selection criteria at different wavelengths 
are able to provide a complete census of such sources. 
Although an unbiased source selection is hard to achieve, 
current surveys in e.g. the obscured Active Galactic Nu- 
clei (AGN) research field have allowed us to assess how 
biases work in source selection at different wavelengths, 
e.g. in the mid-IR vs. X-rays (e.g., Fiore et al. 2008, 



2009; Donley et al. 2008 and references therein). Sim- 
ilar conclusions can be easily extended to other source 
classes and bands. 

If, on the one hand, this approach is clearly consuming 
in terms of telescope/facility time and man-power, on the 
other hand it allows for a proper characterization of the 
intrinsic properties of a given source population. 

Although the science goals which motivate surveys at 
different wavelengths may be different, there are some 
common aspects which will be discussed in the follow- 
ing in order to understand the complex processes behind 
the build-up of source catalogs and the achievement of 
the published results. 



2. SOURCE ASSOCIATION AND IDENTIFICA- 
TION: EXAMPLES FROM THE X-RAY SUR- 
VEYS IN THE COSMOS FIELD 



Motivated by the experience gained over the last decade 
in the field of X-ray surveys by members of our collab- 
oration, hereafter some practical examples and problems 
related to the process of X-ray source identification will 
be presented. In particular, we will mention some issues 
concerning the recent surveys carried out with XMM- 
Newton and Chandra in the Cosmic Evolution Survey 
(COSMOS) field, keeping in mind that some problems 
can be generalized to other surveys and wavelengths. 

The COSMOS surve>Q (Scoville et al. 2007a), an HST 
Treasury Project, is a deep and wide extra-galactic survey 
over a contiguous w 2 deg 2 field designed to have suffi- 
cient area to limit the effects of cosmic variance in the 
investigation of e.g. large-scale structures (LSS; Scov- 
ille et al. 2007b; Gilli et al. 2009), and depth to study 
the high-redshift (zwl-3) galaxy and AGN populations. 



http://cosmos.astro.caltech.edu. 
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Among the wealth of multi-wavelength data (not all cov- 
ering the entire area; see the review talk by McCracken, 
these proceedings) in this region of the sky, the X-ray ob- 
servations carried out with XMM-JVwton (1.55 Ms) and 
Chandra (1.8 Ms) - the so-called XMM-COSMOS and 
C-COSMOS surveys (Hasinger et al. 2007; Elvis et al. 
2009) - cover the entire 2 deg 2 of the COSMOS field 
and the innermost 0.9 deg 2 , respectively. While XMM- 
Newton observations are able to achieve a 0.5-2 keV 
(2-10 keV) flux limit over 90% of the area of w 1.7 x 

1Q-15 cr g cm -2 g-l (~ 9 3 x cr g cm -2 g-1. 

Cappelluti et al. 2009), making X-ray spectral anal- 
ysis for a significant fraction of sources feasible (e.g., 
Mainieri et al. 2007), Chandra allows for a factor « 10 
deeper X-ray flux detections in both the soft and hard 
bands. These depths are well suited for both AGN 
and normal/starburst galaxy studies in the X-ray band 
(Ranalli et al., in preparation) and, given the typically low 
background in Chandra exposures, allow detection of 
faint X-ray emission using stacking analysis techniques 
(see Fiore et al. 2009 for a recent application). The fi- 
nal XMM- and C-COSMOS catalogs comprise 1822 and 
1761 point-like X-ray sources above a detection likeli- 
hood of 10 and 10.8, respectively (Cappelluti et al. 2009; 
Elvis et al. 2009). 



2.1. Matching source catalogs 

Despite the ultra-deep exposures in the Chandra deep 
fields (North and South; see Hasinger & Brandt 2005 
for a review, and Luo et al. 2008 for recent updates), 
which are focused on observing the same portion of the 
sky (i.e., centered on the same aimpoint) repetitively, the 
mosaics in most of the current X-ray surveys are based 
on observations of neighbouring regions with large over- 
laps among the different observations (but, for a different 
observing strategy, see AEGIS; Laird et al. 2009). Both 
XMM-COSMOS and C-COSMOS surveys are based on 
a tiling scheme which, in the latter survey, allows us to 
achieve a remarkably uniform exposure over the inner 
0.45 deg 2 (for details, see Elvis et al. 2009 and Puccetti et 
al., in preparation). This translates directly into a sharply 
defined flux limit. In such a complex tiling strategy the 
same source is detected in more than one pointing with 
largely different Point Spread Functions (PSFs), possibly 
causing a worse source position reconstruction, hence a 
larger positional error in the final catalog, despite the po- 
tentially sharp Chandra PSF (FWHM« 0.5") for on-axis 
source detection. 

The first drawback of having a few arcsec "error box" as- 
sociated to each X-ray source, depending on the source 
position across the field-of-view and the observed num- 
ber of counts, is related to the match of the X-ray 
sources to the counterparts observed at other wave- 
lengths. This has strong consequences in the process 
of source identification and classification via optical and 
near-IR spectroscopy, making the already difficult task of 
spectroscopy (because of the faint magnitudes of many 
counterparts) quite challenging, since a larger number of 



sources should be observed and identified. 

Having these issues in mind, the X-ray source catalogs 
in both XMM-COSMOS and C-COSMOS surveys have 
been matched to catalogs at other wavelengths using the 
deep and good-quality I, K, and 3.6 /im data (Brusa et al. 
2007; Brusa et al., in preparation; Civano et al., in prepa- 
ration). In particular, using near-IR data allows for a more 
complete source association, since a fraction of X-ray 
sources, especially at the very faint flux limits probed 
by Chandra, comprises obscured, optically faint (because 
obscured and at high redshift) AGN, as expected from the 
AGN synthesis models of X-ray background (XRB; e.g., 
Gilli, Comastri & Hasinger 2007). 

The adopted source matching method consists of the ex- 
tensively used likelihood ratio (LR) technique (Suther- 
land & Saunders 1992; for some recent applications of 
this method, see Ciliegi et al. 2003; Brusa et al. 2005; 
Laird et al. 2009; Motch, these proceedings). For a given 
optical/near-IR counterpart, LR is defined as the ratio be- 
tween the probability that the source is the correct iden- 
tification and the corresponding probability for a back- 
ground, unrelated object. The chosen threshold for LR 
should therefore be a a careful compromise between com- 
pleteness (i.e., most of the sources are finally recovered) 
and reliability (i.e., the number of possible spurious asso- 
ciations must be kept as low as possible). Without enter- 
ing in too many details about this procedure, we note that 
this method has produced reliable source associations 
for the X-ray sources in both XMM- and C-COSMOS 
surveys. For instance, « 80% of the XMM-COSMOS 
sources had already a counterpart using the I and K-band 
catalogs, and this fraction increased by w 14% using 
Chandra sharper positions. The Spitzer IRAC 3.6 /im 
data increased the number of associations by a further 9% 
and 5% for the XMM- and C-COSMOS sources, respec- 
tively. This means that the final fraction of X-ray sources 
without an optical/near-IR match (i.e., cases where either 
there are multiple sources in the X-ray "error box" above 
the adopted LR threshold or at least one possible counter- 
part below the LR threshold) is 1 1% in XMM-COSMOS 
(Brusa et al., in preparation) and only 1% in C-COSMOS 
(Civano et al., in preparation), thus allowing for statisti- 
cally reliable studies of the multi-wavelength properties 
of the X-ray sources, including those with faint optical 
counterparts. 

Besides the overall good matching results, we note that 
the adoption of the LR technique requires significant ef- 
forts, mostly aimed at verifying the reliability of source 
associations and identifications beyond their "statistical" 
validity. The complexities and subtleties of such proce- 
dure have been faced - at least in the COSMOS survey 
- using visual inspection of the candidate counterparts 
on a source-by-source basis, requiring significant man- 
power on month timescales, given the large number of 
sources available and the continuous updates in the ref- 
erence catalogs, which are extensively used before their 
final releases. 

In this context, Virtual Observatory (VO) can be of pri- 
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mary interest for people working in the field of galac- 
tic and extra-galactic surveys, provided that the avail- 
able multi-wavelength catalogs are as uniform as possi- 
ble (e.g., for photometric -redshift analyses and spectral 
energy distribution studies, magnitudes within the same 
extraction region should be homogeneously selected) and 
rigorously checked. If such conditions were satisfied, 
a further implementation to current VO matching tools 
could be the selection of one or more matching methods 
within VO itself, and running tests of source reliability, 
possibly using simulations to reassure about the obtained 
results. Given the way people are organized in present 
surveys (i.e., each group is mostly involved in reducing 
and analyzing data at a given wavelength, and finally in 
building up a source catalog), it must be said that imple- 
menting source-matching tasks in VO could be of lim- 
ited use. On the other hand, given the large amount of 
data that will be available in the near future using the 
next generation of ground- and space-based facilities, this 
issue will become more and more important and time- 
consuming, thus invoking for a significantly different ap- 
proach in the way we currently deal with observational 
data. 



3. COMING SOON: A QUICK LOOK AT SOME 
OF THE NEW FACILITIES 



In the following, some of the main next-generation facil- 
ities in the survey field will be outlined, along with po- 
tential problems in data handling and future catalog con- 
struction. These facilities will provide powerful tools for 
source identification on relatively short timescales. 



3.1. eROSITA 

Over the next few years, one of the major advances in the 
field of X-ray surveys will be provided by the eROSITA 
mission (Predehl et al. 2007) which is scheduled for 
launch in 20120 This mission is focused on X-ray sur- 
veys, being the natural continuation of the former all-sky 
survey performed in the soft band by the ROSAT satellite 
in early 90' (RASS; Voges et al. 1999). With an estimated 
PSF Half Energy Width (HEW) of « 15" for on-axis ob- 
servations (and average PSF HEW of w 25"over the en- 
tire field-of-view of 0.83 deg 2 ), eROSITA is primarily in- 
tended to perform a 4-yr all sky survey, reaching a flux 
limit of » 1.0 x 10~ 14 (« 1.2 x 10~ 13 ) erg cm" 2 s" 1 
in the 0.5-2 keV (2-10 keV) band, with typical expo- 
sures of « 2 ks, and a 0.5-yr deep survey over 400 deg 2 
at flux limits about a factor w 5 and w 13 fainter in 
the two bands, respectively (with an average exposure of 
«27ks). 

For what concerns the study of the AGN physics and evo- 
lution, the planned all-sky survey will allow for detec- 
tions of w 4 million AGN, 25000 (2100) of which being 



above redshifts z=3 (4). These numbers derive from the 
most up-to-date XRB synthesis models (Gilli et al. 2007) 
and calibrations from the eROSITA team; in particular, 
the on-flight background (cosmic plus particle) has been 
simulated for a L2 orbit. In survey mode, where source 
statistics is photon limited, the sensitivity has been es- 
timated by imposing that a source can be detected if at 
least 9 net counts are present and by assuming a tipical 
power-law spectrum with photon index T=1.8 and Galac- 
tic column density. Then the source counts have been 
converted into a flux using the response matrices avail- 
able at the time of writing. In deep survey mode, the 
background counts have been accumulated in 30" circu- 
lar apertures, and limiting counts have been estimated by 
imposing to have a 3a detection above the background. 
Counts were then converted into fluxes by using the same 
X-ray spectral model as above. 

This approach will represent a major improvement in the 
field of statistical studies of AGN; for comparison, less 
than one thousand AGN are detected in the X-ray ultra- 
deep exposures in the CDF-N and CDF-S, and of the or- 
der of a couple of thousands in COSMOS (see ji|2jT). The 
eROSITA deep survey, although on a limited region, will 
allow for detections of X-ray faint AGN at high redshifts, 
which are currently poorly sampled by current X-ray sur- 
veys (e.g., see Fig. 3 of Vignali et al. 2005). 

The large amount of sources detected by eROSITA will 
require, at some level, multiple campaigns of photo- 
metric and spectroscopic follow-up observations. Just 
to give an idea of the numbers involved in these 
follow-ups, at the conservative 0.5-2 keV flux limit of 
«4x 10~ 15 erg cm -2 s _1 for the eROSITA deep sur- 
vey, we find w 570 sources in XMM-COSMOS, with the 
faintest sources having an I-band (3.6 /im) magnitude of 
24 (22; both magnitudes are in the AB system). Cov- 
ering 400 deg 2 at such depth is currently hard to accom- 
plish, and the whole spectroscopic coverage will be even 
more complex and time-consuming. Given the eROSITA 
PSF, the source identification will not be an easy process 
and will take advantage of follow-up observations, espe- 
cially in the near-IR, where the planned WISE and VISTA 
surveys will provide a significant contribution to source 
identification, hence to observational cosmology. These 
and other planned facilities that will be available to the 
community over the coming years are briefly reviewed in 
the following. 



3.2. WISE 

The Wide-field Infrared Survey Explorer (WISE) is a 
NASA missiorQ (launch scheduled for mid 2009) aimed 
at surveying the entire sky in four bands (3.3, 4.7, 12, 
and 23 /im), with sensitivities of 0.12, 0.16, 0.65 and 
2.6 mJy (after eight repetitions of 8.8s each), respectively. 
Such sensitivities are worse than those achieved by the 
instruments onboard the Spitzer satellite, but are a factoi- 
ds 1000 better than those obtained with the IRAS surveys 



2 http://www.mpe.mpg.de/projects.html#erosita. 



3 http://wise. ssl.berkeley.edu. 
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(Mainzer et al. 2006). Furthermore, although the res- 
olutions allowed by WISE are typically worse (FWHM 
of 6" and 12" in the near-IR and mid-IR, respectively) 
than Spitzer, the large field-of-view of WISE (47 arcmin 2 ) 
will allow to perform a survey of the entire sky in just 
six months. This means that studies like those regard- 
ing obscured high-redshift AGN by means of mid-IR-to- 
optical flux ratio selection, as those carried out recently 
by several authors (e.g., Daddi et al. 2007; Fiore et al. 
2008, 2009; Dey et al. 2008; Lanzuisi et al. 2009), will be 
possible over much larger sky regions (although at higher 
flux limits), provided that coverage at other wavelengths 
(to select and properly characterize the sources) will be 
granted. 

Similarly, detection of obscured AGN from shallow and 
moderate-depth X-ray surveys up to high-redshifts will 
be possible, as shown in Fig. [T] where the optical, in- 
frared and sub-mm spectral energy distribution (SED) 
of a luminous Type 2 quasar at z=1.957 detected by 
the HELLAS2XMM survey (Fiore et al. 2003) is pre- 
sented. The intense star formation derived for this source 
(~ 1500 M Q yr _1 ) is suggestive of an AGN passing 
through a complex phase, being caught presumably be- 
fore the process of expelling the obscuring gas (thus 
quenching the star formation) and becoming a broad-line 
unobscured quasar has started. AGN with similar prop- 
erties do not seem to be uncommon at high redshifts, as 
shown in recent work (e.g., Mainieri et al. 2005; Polletta 
et al. 2008). Through proper broad-band SED fitting, 
WISE and Herschel will be able to constrain the individ- 
ual components (stellar, AGN, and starburst) and estimate 
the nuclear bolometric luminosity (see Vignali et al. 2009 
for details). 



3.3. VISTA 

The Visible and Infrared Survey Telescope for Astron- 
omy (VISTAfl is a 4-m wide-field (1.65 deg diameter) 
survey telescope equipped with a near-IR camera facility 
(at least at the beginning of the operations, although it is 
also capable of operating in the visible with a suitable 
camera), covering the 0.85-2.3 /im wavalength range 
with five broad-band filters (Z, Y, J, H, and K s ) and 
a narrow -band filter centered at 1.18 /im. Located at 
the ESO Cerro Paranal Observatory in Chile, the VISTA 
planned capabilities are such to suitably match the needs 
of the 8 m-class telescopes of VLT. The large field-of- 
view betrays that VISTA is mainly designed for carry- 
ing out large-area surveys in the near-IR domain. In this 
regard, it is properly suited for follow-up identification 
programs (e.g., see § 13.11 ). One of the proposed sur- 
veys is the VISTA Hemisphere Survey that will image 
w 20000 deg 2 in the J and K s bands. About 40% of 
this survey, dedicated to a deep observational program 
searching for the dark energy (Dark Energy Survey - 
DES), will cover sa 5000 deg 2 also in the H band over 
initially 125 nights down to 10cr AB magnitude limits 



T 



0.1 r AGN 




Observed-frame wavelength (/im) 



Figure 1. Observed-frame SED of an X-ray selected, lu- 
minous and obscured quasar from the HELLAS2XMM 
survey (Vignali et al. 2009) at z=1.957. The "total" 
SED (solid line) is the summed contribution of a stel- 
lar component (dotted line), a torus component (long- 
dashed line), and a starburst component (short-dashed 
line). The sensitivity and confusion limits expected for 
WISE (four near-IR and mid-IR bands) and Herschel 
( three mid-IR/far-IR bands with PACS and three far-IR 
bands extending from fa 200 to 670 pm with SPIRE) are 
shown as horizontal bars, whose length is indicative of 
the band covered by each filter or instrument ( as in the 
case of SPIRE). 

in the J, H, and A" s bands of 20.4, 20.0, and 19.4, re- 
spectively (Banerji et al. 2008). Furthermore, the Ultra- 
VISTA survey is aimed at imaging the COSMOS field 
down to unprecedented depths, taking advantage also of 
the deep 3.6^1.5 coverage with the planned Spitzer 
Extended Deep Survey (SEDS) "warm-mission" project. 
The UltraVISTA survey will use all of the near-IR fil- 
ters, including the narrow-band one to search for Lya 
emitters at very high redshifts (z rs 8.8). Finally, we 
mention also the VISTA Deep Extragalactic Observations 
Survey (VIDEO), which will cover 15 deg 2 of the sky 
(some of which already included in SWIRE) with the 
Z, Y, J, H and A" s filters, jointly with the recently ap- 
proved Spitzer Extragalactic Representative Volume Sur- 
vey (SERVS) project. 



3.4. Pan-STARRS 

The Panoramic Survey Telescope and Rapid Response 
System (Pan-STARRS) consists of four telescopes of 
1.8m each located in Hawaii (the prototype PS1 being al- 
ready on site)0 Its innovative design for wide-field imag- 
ing, which comprises a 64x64 CCD array (600x600 
pixel each, for a total of 1.4 billion pixels), is such to 
allow observations of the available sky (w 30000 deg 2 ) 
several times each month, down to r w 24.0 (5a). Al- 
though the immediate goals of Pan-STARRS are related 



http://www.vista.ac.uk/index.html. 



s http://pan-starrs. ifa.hawaii.edu/public. 
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to discover and characterize Earth-approaching objects, 
the huge amount of images guaranteed by this facility 
will allow for fast photometric observations of large sky 
regions (up to as 3000-6000 deg 2 in one night). Further- 
more, the wavelength range of the available filters (g, r, i, 
z, and y, covering the as 4700-10500 A interval) is suited 
to provide a first-order characterization of the sources, 
including an estimate of their redshift via photo-z tech- 
niques. 



3.5. LSST 

The Large Synoptic Survey Telescope (LSST0 repre- 
sents one of the most revolutionary ideas for the next 
decade. Aimed at probing dark energy and dark matter, 
mapping the Solar system and the Milky Way, and explor- 
ing transient object, LSST will offer unique opportunities 
for large-area surveys. It consists of a large, wide-field 
ground-based system designed to obtain multiple images 
covering the as 30000 deg 2 of sky which are visible from 
Cerro Pachon in Chile. The current baseline design is 
characterized by an 8.4m (6.5m effective) primary mir- 
ror with a 9.6 deg 2 field-of-view and a 3.2-billion pixel 
camera (Gilmore et al. 2008), and allows observations 
of about 10000 deg 2 of sky using pairs of 15s exposures 
in two of the six photometric bands (u, g, r, i, z, and y) 
every three nights on average, with typical 5er depth for 
point-like sources of r as 24.7 (Ivezic et al. 2008). The 
"deep-wide-fast" survey mode will be assigned 90% of 
the time, thus allowing for observations of as 20000 deg 2 
in the six bands (over the ten years of operations) and 
reaching a final depth of r as 27.7. The remaining 10% 
of the time will be allocated to special programs such as 
the "Very Deep and Fast time domain survey", reaching 
r as 28. Three billion galaxies are expected to have reli- 
able photometric redshifts at the end of the observations, 
and 20-80 million AGN should be detected (i.e., two or- 
ders of magnitude larger than the current AGN census); 
a fraction of AGN (depending on their absolute magni- 
tudes, redshifts, and observed time intervals) will be rec- 
ognized as such on the basis of their variability, thanks to 
the multiple observing-epoch strategy. 



3.6. Spectroscopic facilities: HETDEX and SDSS- 
III 

While the previous facilities will focus on imaging ob- 
servations of large sky regions, here we present a cou- 
ple of new facilities which will be available to the com- 
munity to carry out spectroscopic programs. The first 
is the HETDEX (Hobby-Eberly Telescope Dark Energy 
Experiment)/VIRUS (Visible Integrated-Field Replicable 
Unit Spectrographs) projectjj which consists on an inno- 
vative set of spectrographs (for a total number of as 34000 
fibers) to observe thousands of galaxies each night. The 

6 http://www.lsst.org/lsst. 
7 http://hetdex.org/hetdex/virus.html. 



HETDEX/VIRUS survey, with the primary goal of ex- 
ploring the issues related to the dark energy through a 
three-dimensional map of a large volume of space, is 
scheduled for 2010-2013 and is going to provide red- 
shifts for about 1 million galaxies. 

Over the 2009-2014 period, the Sloan Digital Sky Sur- 
vey III (SDSS-Hlfl will carry out a program of four sur- 
veys on three scientific themes: dark energy, dynamics 
and chemical evolution of the Milky Way, and studies of 
planetary systems. In particular, the Baryon Oscillation 
Spectroscopic Survey (BOSS) project, in order to map the 
spatial distribution of luminous galaxies and quasars, will 
survey 10000 deg 2 and provide redshifts for as 1 .5 million 
galaxies up to z=0.7, and 160000 quasars at zas2.2-3, us- 
ing a 1000-fiber spectrograph with resolution R«2000. 



4. THE POSSIBLE ROLE OF THE VIRTUAL OB- 
SERVATORY 



Given the large amount of data that will be available to 
the community in the future, it seems plausible that VO 
can have an effective role in data handling and related 
tasks, as actually presented by many authors during this 
workshop. As pointed out in i j2.ll large multi-wavelength 
datasets require homogeneous data selection and valida- 
tion, and quality checks. Having said this, it is imperative 
for all of the people releasing reduced data and corre- 
sponding source catalogs to carry out high-level quality 
checks and exhaustive accompanying information. 

From an X-ray perspective, we have shown that the pro- 
cess of building up and cross-correlating catalogs using 
proper matching criteria (e.g., the likelihood ratio tech- 
nique) is time-consuming and requires, at some level, 
visual inspection to support the statistical reliability of 
the matched sources. If, on the one hand, some checks 
are clearly beyond the possibilities of VO, on the other 
hand data homogenization and automatization of some 
matching procedures (and corresponding checks) would 
be highly appreciated by the community. Such imple- 
mentations would allow for a quicker, more extensive and 
exhaustive exploitation of the available datasets. Sim- 
ilarly, automated tasks to perform stacking analyses at 
different wavelengths and estimate photometric redshifts 
using broad-band matched catalogs would greatly en- 
hance the use of VO, especially by people not directly (or 
deeply) involved in all the issues and problems related to 
surveys and data handling. 

Synergy with VO is undoubtedly a challenging issue over 
the coming years when billion-object catalogs will be 
available and the necessity of optimizing the many-sided 
scientific outputs on relatively short timescales will be- 
come more and more pressing. 



8 http://www.sdss3.org. 
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